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ABSTRACT 

We present new spectroscopic observations for a sample of C(N)-type red 
giants. These objects belong to the class of Asymptotic Giant Branch stars, ex- 
periencing thermal instabilities in the He-burning shell (thermal pulses) . Mixing 
episodes called third dredge-up enrich the photosphere with newly synthesized 
^^C in the He-rich zone, and this is the source of the high observed ratio between 
carbon and oxygen (C/0 > 1 by number). Our spectroscopic abundance esti- 
mates confirm that, in agreement with the general understanding of the late evo- 
lutionary stages of low and intermediate mass stars, carbon enrichment is accom- 
panied by the appearance of s-process elements in the photosphere. We discuss 
the details of the observations and of the derived abundances, focusing in partic- 
ular on rubidium, a neutron-density sensitive element, and on the s-elements Sr, 
Y and Zr belonging to the first s-peak. The critical reaction branching at ^^Kr, 
which determines the relative enrichment of the studied species, is discussed. 
Subsequently, we compare our data with recent models for s-processing in Ther- 
mally Pulsing Asymptotic Giant Branch stars, at metallicities relevant for our 
sample. A remarkable agreement between model predictions and observations is 
found. Thanks to the different neutron density prevailing in low and intermediate 
mass stars, comparison with the models allows us to conclude that most C(N) 
stars are of low mass (M < 3 Mq). We also analyze the ^^C/^^C ratios measured, 
showing that most of them cannot be explained by canonical stellar models. We 
discuss how this fact would require the operation of an ad hoc additional mixing, 
currently called Cool Bottom Process, operating only in low mass stars during 
the first ascent of the red giant branch and, perhaps, also during the asymptotic 
giant branch. 
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Subject headings: nucleosynthesis — stars: abundances — stars: AGB — stars: 
carbon 

1. Introduction 

The production of elements heavier than iron is normally ascribed to neutron captures. 
At the low neutron density prevailing in Asymptotic Giant Branch (AGB) stars, this 
nucleosynthesis process is referred to as the s-process, where the neutron-capture path 
proceeds along the valley of P stability and most unstable isotopes encountered along the 
s-path preferentially /?-decay instead of capturing a further neutron. However, in a number 
of cases, competition between neutron capture and beta decay leads to a branching in the 
synthesis path, which is determined by the neutron density and the neutron capture cross 
section of the unstable nucleus and its /9-decay half-life (see e.g. Kappeler, Beer, & Wisshak 
1989). 

The discovery of the unstable element technetium in the spectra of AGB stars by 
Merrill (1952) clearly showed that the s-process actually occurs in these red giants of 
low to intermediate mass, 1.5 < M/Mq < 8. The particular structure of AGB stars (two 
alternate burning shells, one of H and an inner one of He, surrounding an inert degenerate 
CO core) fulfills the necessary physical conditions for the s-process nucleosynthesis to 
occur. It is in the He-rich zone between the two burning shells (hereafter He intershell) 
that s-elements are produced by slow neutron captures on seed nuclei. Periodically the He 
intershell is swept by a convective instability induced by a He-burning runaway (thermal 
pulse, TP) where ^^C is synthesized. After each thermal instabihty, newly synthesized ^^C 
and s-elements are mixed with the envelope by a mixing episode called third dredge-up 
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(TDU) (Iben & Renzini 1983). The repeated action of TPs and TDU episodes during the 
AGB phase makes the star to eventually become carbon-rich, i.e. showing a ratio C/0 > 1 
in the envelope. Abundance studies of C-stars (N and SC types) in the Galaxy have indeed 
confirmed this figure (Utsumi 1985; Abia & Wallerstein 1998). 

One debated subject concerning the s-process nucleosynthesis in AGB stars is the 
knowledge of the neutron source. Two competing sources have been envisaged: the 
^^C(Q;,n)^^0 and the ^^Ne(Q;, n)^^Mg reactions. During the last decade a number of studies 
in AGB stars of different spectral types have shown that the s-element abundance pattern 
can only be understood if the ^^C(q;, n)^^0 reaction provides the bulk of the neutron flux 
at low neutron densities (A^„ < 10^ cm~^) and relatively low temperatures (T ~ 1 x 10^ K) 
(cf. Wallerstein et al. 1997; Busso, Gallino, & Wasserburg 1999, and references therein). In 
AGB stars, this reaction occurs in radiative conditions during the interpulse phase, in the 
top layers of the He intershell (Straniero et al. 1995, 1997; Herwig et al. 1997; Gallino 
et al. 1998; Herwig 2000). In fact, the convective instability powered by a TP spreads 
over the whole He intershell the ^^C produced by the 3a reactions. When TDU takes place, 
a few amount of protons may diffuse from the H-rich envelope down to the He intershell. 
Then, during the interpulse, the region where these protons are stored heat up so that ^^C is 
rapidly formed, via proton captures on ^^C, and later on, when the temperature approaches 
10® K, the ^^C(Q!,n)^^0 neutron source is activated (Busso, Gallino, & Wasserburg 1999). 

On the other hand, the ^^Ne(Q;, n)^^Mg neutron source requires higher temperatures 
(T > 3.5 X 10® K), as those currently found in the TP of intermediate mass models (i.e. M 
> 4 Mq) (see e.g. Straniero et al. 2001). The maximum temperature achieved in low mass 
AGB stars (LMS) (M < 3 M©) instead does not exceed T ~ 3 x 10® K (Iben & Truran 
1978; Straniero et al. 1997). Therefore, in LMS the ^^Ne neutron source is only marginally 
activated, and the ^^C source is the main supplier of neutrons. Conversely, in AGB stars of 
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M > 4 Mq, both the ^^Ne and the ^^C neutron sources are efficiently activated, during the 
thermal pulse and during the interpulse, respectively. 

The neutron burst by the ^^Ne neutron source occurs typically at high peak neutron 
densities, Nn > 10^^ cm~^ (see Vaglio et al. 1999; Busso, Gallino, & Wasserburg 1999). 
This strongly favors the production of the neutron-rich nuclides affected by the branching 
in the s-path, like ^^Kr, ^^Rb and ^^Zr. Therefore, a different s-element pattern is expected 
depending on whether one or the other neutron source is more active, thus allowing us to 
estimate the neutron density at the s-process site and consequently the initial AGB mass. 
The typical mass of the C-stars may be inferred by means of their luminosity distribution. 
In the LMC the C-stars luminosity function is peaked at Mboi ~ —4.5 (see e.g. Hughes & 
Wood 1990). A higher value (Mboi ~ —5) is found in the Milky Way. These luminosities 
are typical of low mass AGB stars. Very few C-stars have been found at Mboi < —6 (van 
Loon et al. 1998, 1999; Wallerstein & Knapp 1998; Trams et al. 1999), the mass of C-stars 
in the Magellanic Clouds being estimated to be lower than 2.7 Mq (Progel, Mould, & 
Blanco 1990). Such an occurrence seems to exclude that stars with M > 4 Mq could 
become C-stars. However, it may be possible that the lack of bright C-stars is due to an 
observational selection effect: the brightest AGB stars could be obscured by an optically 
thick wind before they reach the C/0 > 1 phase. On theoretical grounds, the question of 
the C-star mass has been a debated subject, because of the difficulty of finding TDU in 
LMS (Iben 1981). This problem has by now been overcome, since the TDU occurrence has 
been found in AGB stars with initial mass as low as 1.5 Mq (Lattanzio 1989; Straniero et 
al. 1997; Herwig 2000). 

In addition, AGB stars of intermediate mass develop a large core mass (Mh > 0.9 Mq), 
and the temperature attained at the base of the convective envelope during the interpulse 
period becomes larger than 70 — 80 x 10^ K, so that ^^C is converted into ^''N. This is the 
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so-called Hot Bottom Burning (HBB, Sugimoto 1971; Iben 1975; cf. Iben & Renzini 1983) 
which is currently found in the most recent computations of intermediate mass AGB models 
(e.g. Blocker & Schonberner 1991; Boothroyd, Sackmann, & Ahern 1993; Lattanzio et al. 
1996; Forestini & Charbonnel 1997; Straniero et al. 2001). Such an occurrence could explain 
the lack of bright C-stars. However, the HBB efficiency also depends on the envelope mass. 
When, as a consequence of the mass loss, the mass of the envelope is reduced to ~1.5 Mq, 
HBB eventually ceases, while TDU continues to mix carbon from the He intershell with the 
surface (Boothroyd & Sackmann 1992; Frost et al. 1998). This could lead to the formation 
of a bright, possibly obscured, C-star. Note that during the HBB phase the ^^C/^^C ratio 
would remain very low and close to its equilibrium value of 3.5, and eventually rise up to 
5 — 10 when C/0 ~ 1 is achieved in the envelope (Lattanzio & Forestini 1998). 

In this study abundances of Rb, Sr, Y and Zr in a sample of Galactic C-stars are 
simultaneously derived for the first time. By comparing these findings with those predicted 
by s-process nucleosynthesis models at the ^^Kr-branching in AGBs of different masses, 
we infer about the stellar mass of C-stars and the neutron source playing the main role 
in the s-process. In § 2 we describe the properties of the ^^Kr-branching. § 3 presents 
the observations and analysis. In § 4 the abundances derived in our sample of C-stars are 
discussed and compared with the s-elements abundance pattern found in others AGB-stars. 
Finally, in § 5 and § 6 we discuss the observational results in the framework of the current 
s-process nucleosynthesis models in LMS and IMS. Our conclusions are summarized in § 7. 

2. The reaction branching at ^^Kr and the neutron density 

The peculiarities of the neutron fiow near the magic neutron number N — 50, and 
through the unstable ^^Kr in particular, offer the possibility of making an observational 
estimate of the neutron density. Neutron capture on ^''Kr leads either to the 10.7 yr ground 



- 7- 



state of ^^Kr or to the 4.5 h isomeric state at 350 keV, which /3-decays roughly by 80% to 
^^Rb and by 20% to the ^^Kr ground state through a 7-transition. The production rate 
of S5j^].iso jg qnite important, amounting to about 50% of the total neutron capture cross 
section on ^^Kr. The half-life of the ground state of ^^Kr is long enough for the s-process 
flow to proceed further to ^^Kr and then to ^^Rb, in competition with /3-decay to ^^Rb, 
depending on the neutron density. When the neutron density is higher than a few 10^ 
cm~^, the neutron channel is open and the s-process path feeds the neutron magic nuclei 
^^Kr and ^^Rb. Since ^^Rb has a neutron capture cross section larger by a factor ~ 10 than 
^^Rb, the total abundance of Rb as well as its isotopic mix is very sensitive to the neutron 
density. Indeed, the total Rb abundance can differ by one order of magnitude depending on 
whether the low or the high neutron density routing is at work (see Beer & Macklin 1989 
for a detailed description of this branching) . Therefore, the relative abundance of Rb to 
other elements in this region of the s-process path, such as Sr, Y and Zr, can be used to 
estimate the average neutron density of the s-process and, as a consequence, to infer the 
mass of the star. In fact, in the last 15 years this kind of study for Rb has been successfully 
performed in field stars (Tomkin & Lambert 1999) and in AGB stars of different spectral 
types: Ba stars (Tomkin & Lambert 1983; Malaney & Lambert 1988), MS and S giants 
(Lambert et al. 1995), yellow symbiotics (Smith et al. 1997) and some s-process- rich stars 
in iij Cen (Vanture, Wallerstein, & Brown 1994; Smith et al. 2000). These abundance 
studies have added further evidence favoring the ^•^C(a,n)^^0 reaction as the neutron 
source. The observed sources in the above studies often belong to binary systems in 
which the companion is now a white dwarf (Jorissen et al. 1993). Thus, they owe the 
s-element overabundances to mass transfer from an ancient AGB star (the progenitor of 
the white dwarf companion) and not to an ongoing s-process. As to C(N) stars, the only 
spectroscopic estimates of s-process elements available so far are those by Utsumi (1985), 
which are however based on low-dispersion photografic plates, nor Rb was measured at the 
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time. With the improved set of spectroscopic data presented here we are able to determine 
the Rb abundance and to provide a much better estimate of Sr, Y, Zr, thus allowing us to 
examine the general problem of the initial mass also for the class of the C(N) stars. 

3. Observations and analysis 

The spectra of the C-stars analyzed in this study were taken at the Calar Alto 
Observatory using the 2.2 m telescope, and at the La Palma Observatory with the 4.2 m 
WiUiam Herschel telescope (WHT). They form part of a more extended study concerning 
the chemical composition of Galactic C-stars: see Abia & Wallerstein (1998), hereafter 
Paper I, and Abia & Isern (2000), hereafter Paper II. The echelle spectra obtained have a 
continuous coverage from A ~ 4000 to 9000 A and a resolution A/ A A = 40000 for the Calar 
Alto spectra, or 60000 for those taken at WHT. The spectra are of good quality with S/N 
ratios from ~ 50 in the blue part of the spectrum, up to more than 100 in the red part. 
The spectra were reduced following standard techniques using the IRAF software package. 
To remove telluric absorptions, the reduced and calibrated spectra of the target stars were 
obtained by subtracting the spectrum of a bright, hot and rapidly rotating star. 

The stellar parameters of the stars studied were derived in a similar way as in Paper II. 
In short: effective temperatures were mainly deduced from the calibration of the observed 
bolometric flux to the infrared flux in the L'— band (3.7/xm) versus T^s (stars marked with 
an asterisk in Table 1), or from the (J — L')^ vs. T^s cahbration by Ohnaka & Tsuji (1996). 
Infrared photometry was taken from Fouque et al. (1992) and Noguchi et al. (1981). When 
photometry was not available we adopted the T^s values derived from angular diameter 
measurements (Dyck, van Belle, & Benson 1996). Otherwise, we adopted the value given 
in Lambert et al. (1986) for some stars in common with this work. We have used the 
grid of model atmospheres for C-stars computed by the Uppsala group (see Eriksson et 
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al. 1984, for details). These models cover the Teff = 2500 - 3100 K, C/0 = 1.0 - 1.35 
ranges and all have the same gravity, log g = 0.0. Considering the T^g of the program stars 
and a typical luminosity 10** Lq, this is a common value for AGB stars. Nevertheless, our 
gravity might be in error by as much as ± 1 dex due mainly to the uncertain estimate of 
the luminosity (Mboi). The CNO abundances in the model atmosphere of a given star were 
taken from the literature when available (Lambert et al. 1986), or estimated from spectral 
synthesis in the wavelength ranges analyzed here. The absolute abundances of CNO that 
can be estimated in that way are rather uncertain; however, the C/0 ratio that mainly 
determines the global shape of the spectrum, is derived with good accuracy (~ ±0.05). 
Whenever possible, *^C/*"^C ratios were obtained from our echelle spectra by synthetic 
fits to ^^CN features in the 7990 - 8040 A range following the same procedure as in Abia 
& Isern (1996). Otherwise, they were taken from the literature (Lambert et al. 1986; 
Ohnaka & Tsuji 1996). Some of the *^C/*^C ratios derived here differ from those in Abia 
& Isern (1997) because of the different T^s values and an updated CN and C2 line list in 
the 7990 - 8040 A range used here. Finally, a typical microturbulence parameter of ^ ~ 2.2 
km s~^ was adopted in the analysis of all sample stars (Lambert et al. 1986). Table 1 
shows the atmosphere parameters for the program stars. 

The metallicity [M/H]* of our sample sources was derived from selected metallic lines 
in the stellar spectra (see Table 3). Most of these lines are located in the wavelength ranges 
4450 — 4650 A and 4700 — 4950 A. In these wavelength windows, molecular absorptions are 
weak in C-stars and metallic lines can be identified more easily. The metallic lines were 
selected according to the same careful criteria described in Paper 1. Most of the /-values 
were taken from the literature or from the VALD (Piskunov et al. 1995) and Kurucz's 

*We adopt here the usual notation [X] = log (X.) program-star - log (X.) comparison-star for 
the abundance of any element X. 
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CD-ROM No. 23 data bases. In some cases, we derived solar /-values using the Holweger 
& Miiller (1974) atmosphere model for the Sun. The abundances of metals were derived 
by the usual method of equivalent width measurements and curves of growth calculated in 
LTE. Upper limits to the equivalent widths (see Table 3) were not considered in deriving 
abundances. The [M/H] value for each star in Table 2 is the mean value obtained from the 
selected metallic lines. 

The abundances of the low-mass s-elements around the ^^Kr-branching (Rb, Sr, Y and 
Zr) were derived from spectroscopic features selected following the same criteria as those for 
metallic lines. Despite the fact that our C-stars show many heavy-element absorption lines, 
only a few of them were considered useful for abundance analysis. The reason of this is the 
severe blending with CN and C2 lines, which is evident even at the high spectral resolution 
used here. The lines selected are shown in Table 3. Spectroscopic parameters were taken 
from the hteraure (see Table 3) since due to their weakness in the solar spectrum, the 
derivation of solar /-values is very uncertain. Nevertheless, our ^'/-values are in very 
good agreement with the solar g( /-values list derived by Thevenin (1989, 1990). For a 
sample of 30 lines in common we found a mean difference of 0.12±0.20 dex. In a few 
cases, however, there is a significant difference which, in our opinion, clearly shows that 
the derivation of solar (//-values from very weak hues is not always safe. For instance, our 
main Y abundance indicator is the Y I line at 4819 A. This feature is quoted with a 3 mA 
absorption in the Solar Spectrum Atlas by Moore, Minnaert & Hootgast (1966). However, 
Hannaford et al. (1982) note that the uncertainty in the wavelength of this hne could be 
as high as 25 mA and very probably, Y being not the main contributor to this absorption 
in the Sun^. In fact, using a solar gf -value for this line (~ 1 dex higher than the value in 

^Indeed, the identification of this line as Y in the Solar Spectrum Atlas of Moore, Minnaert 
& Hootgast (1966) is quoted with a question mark. 
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Table 3) in the analysis of WZ Cas (our reference star, see below), we obtain a very low Y 
abundance ([Y/H]~ —1) which is unrealistic given the near solar metallicity of this star and 
its level of heavy element enhancement (see Paper II) . The same figure is found with two 
additional discrepant Ti I lines with the Thevenin's list (at 4812 and 4821 A): unreahstic 
Ti abundances are obtained in WZ Cas if solar gf -values are used (the absorptions in the 
Sun are 8.5 and 3 mA, respectively). Finally, our gf -value of the Fe I line at 5848 A also 
differ significantly with that by Thevenin. However, the iron abundance derived from this 
line in the only star where it is identified (Z Psc), is in agreement within ±0.2 dex with 
that derived from other Fe lines in this star. 

Abundances of s-elements were also derived from curves of growth in LTE except for 
those fines presenting clear blends, namely: the 7800 A Rb I, 7070 A Sr I and 4815.05, 
4815.65 A Zr I lines, for which we used spectral synthesis (see Table 3). CN and C2 lines in 
these wavelength ranges were kindly provided by B. Plez and P. de Laverny according to the 
spectroscopic predictions by de Laverny & Gustafsson (1998). Final individual abundances 
were obtained as the mean of the abundances derived from each line. Unfortunately, given 
the important error in the measurement of the equivalent widths and the small number 
of lines used, it was impossible to check the atmospheric parameters (Teflf, log g and ^) 
using the requirement that individual abundances derived from lines of different intensity, 
ionization state and excitation energy, have to be nearly equal. 

In order to reduce systematic errors due to wrong gf-values and departures from LTE, 
certainly present in late-type giants (Tomkin & Lambert 1983), we performed a differential 
analysis line-by-fine using WZ Cas (C9,2J; SC7/10) as the comparison star. This star has a 
similar T^e and surface gravity as the stars studied here and its spectrum is also dominated 
by CN and C2 absorptions. However, because the C/0 ratio in its atmosphere is very 
close to unity (C/0 ~ 1.005) the spectrum is less crowded with molecular absorptions. 
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and thus the atomic hne identification is easier. This permits a more accurate abundance 
analysis. Abundance studies of this star reveal that its metallicity is almost solar, [Fe/H] ^ 
0.0 (Lambert et al. 1986; Paper II), with no or very small s-element overabundances: the 
mean heavy element enhancement is [h/Fe]= 0.02 ± 0.20 , where 'h' includes Sr, Y, Zr, Nb, 
Ba, La, Ce, and Pr (see Paper II). All the abundance ratios shown in Table 2 are relative 
to WZ Cas\ 

The derivation of the Rb abundance requires special attention. In C-stars the line used 
here is strongly blended with CN, C2 lines and some atomic absorptions. Therefore, the 
Rb abundance must be derived by spectral synthesis. We have used an updated version 
of the molecular and atomic line hst used in Papers I and II. The new line hst includes 
some additional ^^C^^N, ^^C^^N and ^^C^^C lines computed by B. Plez. We have also 
modified the value of a weak Si I line to the blue of the Rb line according to Tomkin 
& Lambert (1999). The hyperfine structure of the Rb I line was considered adopting a 
Solar System isotopic ratio ^^Rb/^^Rb = 2.59 (Anders & Grevesse 1989). Note that the 
isotopic Rb mixture can certainly be altered as a function of the neutron density in the 
stars. Unfortunately, this isotopic ratio cannot be measured spectroscopically. Nevertheless, 
we checked that other isotopic mixtures will not change significantly our conclusions (see 
below). We have adopted the meteoritic Rb abundance, log (Rb/H) + 12 = 2.40 ± 0.05 
(Anders & Grevesse 1989). We preferred to adopt the meteoritic abundance rather than 

^The [h/Fe] value in WZ Gas differs slightly from that derived in Paper II [h/Fe] = +0.16, 
because of the downward revised Sr, Ba and La abundances here. This new analysis was 
done with a very high resolution spectrum (R~ 180000, S/N> 70) obtained, for another 
scientific purpose, with the 2.5 m NOT and the SOFIN echelle spectrograph at La Palma 
Observatory. In the re-analysis of WZ Gas we have used Sr, Ba and La lines mainly placed 
in the 4800 - 4950 A spectral range. 
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the photospheric one (2.60 ± 0.15) as the reference value because of the weakness of the 
7800 A Rb I hne in the Sun (< 1 mA) from which the photospheric abundance is derived. 
Instrumental and macroturbulence broadening were included in the Rb synthesis. The 
macroturbulence broadening was set by matching the profile of an almost clean nearby Ni I 
hne at 7788.9 A. 

A striking result immediately appears in the analysis of the 7800 A Rb I line: the Rb 
abundances derived are unrealistically low considering the near-solar metallicity of most 
C-stars analyzed here. In fact, the derived abundances of Rb are lower than solar by a factor 
ranging typically between four and ten. This figure has also been found in similar studies on 
AGB stars (Plez, Smith, & Lambert 1993; Lambert et al. 1995; Abia & Isern 2000). The 
reason for this is not well understood, although N-LTE effects due to overionization might 
partially explain this finding. In fact, due to its low ionization energy (4.16 eV), rubidium 
is mostly ionized at the depth of the 7800 A hne formation. Abia, Pavlenko, & de Laverny 
(1999) showed that N-LTE effects in C stars due to overionization in the formation of the 
resonance Li I at 6708 A (another alkali clement with an ionization energy similar to Rb, 
5.39 eV) might extend up to -1-0.6 dex (in the sense of N-LTE minus LTE abundances). 
Surprisingly, the same effect is found when deriving K abundances (ionization energy 4.34 
eV) from a nearby resonance K I line at 7698.70 A; the K abundances derived are also 
unrealistically low: typically [K/H] ~ —1^. No nuclear process able to deplete K in stars is 
known. Although both the Li and K hues are usually stronger than the Rb I hne in C-stars 
and, in consequence, they should form higher in the atmosphere where larger departures 
from LTE are expected, upward corrections of the Rb abundance are not excluded until 

Certainly N-LTE effects are important in the formation of the resonance 7698 A K I line 
in the Sun (Takeda et al. 1996). However, as far as we know, there is no a similar study in 
C-rich giants. 
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a quantitative analysis of N-LTE effects in the resonance Rb line is made. However, we 
believe that N-LTE effects are not probably the full story. In our opinion there also seems to 
be a problem related to the continuum opacity in the 7800 A region. Indeed, the metallicity 
derived from (almost unblended) Ni I and Fe I hues near the Rb I hne is significantly lower 
(by a few tenths of dex) than the one derived from the metallic lines in the blue part of the 
spectrum. This suggests that the continuum absorption coefficient of the real atmosphere 
in the Rb spectral region is not well reproduced by our model atmospheres. Tests made by 
modifying the model atmosphere parameters in a quantity similar to the expected errors 
(see Table 4) do not solve the problem. However, when trying to artificially simulate this 
missing opacity by increasing the H~ opacity (namely, increasing the electronic pressure in 
the model atmosphere), the derived metallic and Rb abundances increase by tenths of dex. 
Again, we have to wait for improved model atmospheres of C-stars (Plez 1999) to study 
this problem better. Therefore, because of the uncertain Rb abundances, we computed the 
[Rb/M] ratio using the metallicity ([M/H]) derived from the metallic lines close to the Rb 
line. The [Rb/M] ratios in Table 2 were computed in that way, also using WZ Cas as the 
comparison star. 

Unlike Rb, strontium has many spectroscopically accessible lines. However, in C-stars 
most of them are very strong and/or located in crowded regions of atomic and molecular 
absorptions. As the main indicator of the Sr abundance we have used the Sr I line at 
7070.1 A. In some of the stars observed at a very high resolution {R ~ 60000), it was 
possible to derive Sr abundances from two additional Sr I fines at 4811 and 4872 A. In 
the neighbourhood of these lines, however, the continuum position is uncertain. The Sr 
I line at 7070 A is blended with some ^^CN lines of moderate intensity, so that spectral 
synthesis is required. We obtain a solar gi/-value for this line from the solar spectrum 
using the photospheric Sr abundance 2.90 ± 0.06 (Anders & Grevesse 1989). No hyperfine 
structure was considered for the 7070 A Sr I line, since this absorption is rather weak 
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(log Wx/X < —4.5) even in C-stars with large Sr enhancements. Due to the blend with 
^^CN hnes, the abundance of Sr derived from this feature is sensitive to the ^'^C/^^C 
ratio adopted. Unfortunately, given the current uncertainties in the derivation of CNO 
abundances and the atmospheric parameters of C-stars, the typical error in the derivation 
of the ^^C/^^C ratio is still high, ±(6 - 13): the higher the ^^C/^^C ratio, the higher the 
error (see Abia & Isern 1996). This introduces an extra uncertainty into the derivation 
of the Sr abundance from this line, mainly in stars with ^^C/^^C < 30. Therefore, the 
[Sr/M] ratios derived here have to be considered with some caution^. On the other hand, 
synthetic fits in the Sr spectral range include a clean Ti I hue at 7065.50 A that was used 
as an additional check to the [M/H] ratio derived from metallic lines in the blue part of the 
spectrum and to set the macroturbulence broadening. 



3.1. Errors 

Errors involved in the derivation of abundances in C-stars are large. The sensitivity 
of the derived abundance ratios to the model atmospheres was assessed by changing the 
atmospheric parameters over a plausible range of T^s, ^, CNO abundances, C/0 and 
i2q^i3q ratios. Obviously, for a given element the change in the abundance created by 
varying each atmospheric quantity depends on the excitation energy and on the intensity of 
the specific line (we have used only neutral lines). In Table 4 we show the sensitivity of the 
[X/H] ratios to the above sources of error. We have also included the error introduced by 
an uncertainty of 5% in the location of the continuum in the spectrum (if spectral synthesis 

^Nevertheless, in the stars where it was possible to derive Sr abundances from the hnes 
in the blue, the dispersion never exceeded ±0.3 dex, which is consistent with the expected 
total error. 
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is used) or due to errors in the equivalent width measurements (when curves of growth 
are used). We do not include the effect of uncertainties in the gravity since we have only 
atmosphere models computed with log g = 0.0. However, since we have used only neutral 
lines in the analysis, our abundances would not be very sensitive to changes in gravity 
within a variation of ±0.5 dex in log g. The values reported in Table 4 were computed using 
equivalent widths and model atmosphere parameters that are representative in the program 
stars, namely: -log W^/X = 4.3 to 5.0, T^s = 2850 K, log g = 0.0, ^ = 2.2 km s-\ C/0 = 
1.05 and ^^C/^^C = 40. Prom this table it is apparent that the mean sources of error come 
from uncertainties in T^s, microturbulence and errors in the equivalent width measurement 
or the spectral synthesis fit (by eye). Indeed, the microturbulence in variable stars might 
differ with wavelength in the spectrum (due to differences in depth of line formation) , and 
might change with the pulsation phase of the star. We believe, however, that an uncertainty 
in the microturbulence as large as ±1 kms^^ can be discarded because this can be checked 
from synthetic fits to CN lines in the spectral regions studied. For instance, synthetic 
fits to the strong and crowded CN lines in the 7070 A region show that such a variation 
respect to our adopted value (2.2 kms~^) would either result in CN hues that are too 
narrow and free of significant blending with adjacent lines, or produce extreme saturation 
effects. In the former case, the convolution parameters required to fit the spectrum would 
be very large. In the later, the resulting spectrum would be very smooth with an unrealistic 
continuum position. Furthermore, the effect of variations in the microturbulence may also 
be distinguished from that of variations in the C/0 ratio since the former deepen and 
become broader and more blended with the neighboring features. We note that Lambert et 
al. (1986), analyzing weak CO and CN lines (log (Wx/X) < —4.5) in a sample of 30 C-stars 
(many of them studied here), obtain a mean value ^ = 2.16 ± 0.20 kms~^. On the other 
hand, note from Table 4 the small sensitivity of the Y and Zr abundances and metallicity 
to uncertainties in the C/0, ^^C/^^C and CNO/H ratios. This is because these abundances 



-17- 



are derived using curves of growth from unblended (with CN or C2 features, as far as our 
molecular hne lists indicate) lines. 

The quoted accuracy of our gf-values range between 5% — 25%, which is translated 
to an error in the derived abundances between 0.05-0.25 dex depending on the intensity of 
the line. However, as all the abundances derived here are set in ratio to WZ Cas, errors 
in the gf -values are not the relevant ingredient. More important is the random error in 
the abundance of elements represented by a few number of lines or just one. In this case, 
the quality of the lines in question is the critical point. For instance, wc consider the 
7800 A Rb I line as a good indicator since our synthetic spectra reproduce quite well this 
spectral region. The same apphes for the Zr I blend at 4815 A and for the Y I hne at 4819 
A, which appears safely clean of blends in the majority of the spectra. In any case, an 
estimate of this random error can be evaluated from the resulting dispersion around the 
mean abundance value in the cases where more than three lines are used. This can be done 
for a few stars (see Table 3) and only for Ti and Fe abundances. We found a dispersion 
in the range 0.1-0.25 dex. On the other hand, systematic errors, not included in Table 
4, may be present affecting the [X/M] ratios. The leading contributor to a systematic 
error is probably our assumption that N-LTE effects are reduced when the program stars 
are analyzed with respect to WZ Cas. Note also that our stars are variable, therefore, 
the effective temperature might change with the pulsational phase. The quoted AT^s in 
Table 4, thus, refers to the typical expected uncertainty in the derivation of this parameter 
according to the different methods used: photometry, infrared flux or angular diameter 
measurements (see above). Focusing on Rb for instance, a Rb isotopic mix different than 
solar would modify the derived Rb abundance. Nevertheless, variations of the ^^Rb/^^Rb 
by a factor of five implies a maximum variation of only ±0.1 dex in the Rb abundance. 
On the other hand, if the solar photospheric Rb abundance is preferred as a reference, the 
[Rb/M] ratios in Table 2 have to be reduced by 0.2 dex. Other sources of error such as the 
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uncertainty in the model atmospheres (hydrostatic piano-parallel approximation, velocity 
stratifications, poly-atomic opacities), dust or chromospheres are at present impossible to 
quantify. We estimate a total typical error (non systematic) of ±0.30 dex for [M/H], ±0.30 
dex for [Rb/H], ±0.35 dex for [Sr/H], ±0.30 dex for [Y/H], and ±0.40 dex for [Zr/H]. 
Adding-up quadratically to these figures the random error due to the use of a few number of 
lines, the total error in the [X/H] ratios would range between ±0.3-0.45 dex. Errors in the 
[X/M] and [X/Y] ratios are, in general, lower since these ratios are less sensitive than the 
[X/H] ratios to the atmosphere parameters: sources of error affecting in the [X/H] ratios to 
the same sense cancel out when computing these element abundance ratios. 

Below, we compare our results with previous abundance analyses and with theoretical 
predictions of s-process nucleosynthesis in C-stars of different masses. 



4. The abundance ratios and their characteristics 

A detailed chemical analysis of heavy element abundances from Rb to Ce in a larger 
sample of Galactic C(N)-stars will be presented in a forthcoming paper (Abia et al. 2001, 
in preparation). Here we restrict the discussion to the low-mass s-elements (Rb, Sr, Y 
and Zr). Table 2 shows the abundance ratios with respect to the metaUicity derived in 
our stars referred to WZ Gas. Most stars are of near-solar metaUicity, except V GrB and 
lY Hya, which can be considered as metal-poor G-stars (see below). The stars present 
moderate light s-element (Zr, Y, Sr) enhancements, namely: <[Rb/M]> = 0.27 ± 0.24, 
<[Sr/M]> = 0.41 ± 0.19, <[Y/M]> = 0.59 ± 0.17 and <[Zr/M]> = 0.54 ± 0.22. These 
overabundances are significantly lower than those reported by Utsumi (1985) in a similar 
analysis. Excluding Rb (not studied by Utsumi), our [Is/M] ratios, with 'Is' = (Sr, Y, Zr), 
are 0.4 — 0.6 dex smaller. However, as discussed in Paper H, the abundances derived by 
Utsumi (1985) were based on atomic line identifications made in a previous work (Utsumi 
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1970), and indeed most of the heavy-element hnes identified in this work appear as clear 
blends at the spectral resolution used here (a factor three higher than Utsumi's one). This 
necessarily leads to abundance overestimations. Additionally, Utsumi used in his analysis 
the Allen's (1976) Solar System abundances as reference values, which are lower by 0.05 
(Sr), 0.44 (Y) and 0.1 (Zr) dex than the reference values used here (Anders & Grevesse 
1989). This explains part of the differences found with respect to Utsumi's analysis. 

In any case, the light s-element abundance ratios derived here are lower than previously 
thought and contrast with the higher overabundances found in other C-rich s-enhanced 
stars or their descendants (Ba , CH and SC-stars: see for instance Vanture 1992a,b; Smith 
& Lambert 1988), and in fact are more similar to the findings of 0-rich S stars (Smith & 
Lambert 1990). We remind that during the AGB phase it is commonly assumed that the 
abundance of ^^C and s-elements increases in the envelope of the star along the spectral 
sequence M— >^MS— >S— >SC— >C. As the TP and TDU episodes occur repeatedly, and since 
most of our stars have C/0 ratios near unity, they should show abundances similar to 
those found in SC-stars^. However, comparison with the average low-mass s-element 
enhancements found in SC-stars (see Paper I) shows a difference of about -1-0.5 dex, in the 
sense of SC-stars minus C(N)-stars. We believe that this difference is due to several reasons. 
First, the abundances in Paper I were derived from lower resolution spectra (a factor ~ 2). 
Second, despite the fact that SC stars show less intense molecular absorptions (CN & C2), 
due to their C/0 ratio being very close to unity, the abundance analysis in Paper I was 
based on atomic lines identified in spectral ranges (5900 — 8000 A) that were much more 
crowded with molecular absorptions than those mainly used here (4750 — 4950 A). These two 

^SC-stars are AGB stars with a C/0 ratio very close to unity within 1% or less. Further- 
more, they show spectroscopic characteristics which are used to distinguish them from the 
C(N)-stars (see Keenan & Boeshaar 1980, for a complete definition). 
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facts might have produced systematic abundance overestimates due to undetected blends. 
However, the main reason is probably due to differences in the model atmospheres used. In 
Paper I atmosphere models for SC-stars from D. Alexander (private communication) were 
used. We checked that for the same atmosphere parameters, the thermal structure in the 
higher layers differs significantly from that obtained in the Uppsala models used here. In 
fact, the abundances derived using Alexander's models differ from those derived using the 
Uppsala models by 0.4 to 0.7 dex!, in the sense of Alexander minus Uppsala abundances: the 
lower the excitation energy of the line, the higher the abundance difference is. This affects 
the s-elements abundance estimates more than the metallic abundances, since the latter are 
derived from higher excitation energy lines. If this effect were considered, the [X/M] ratios 
in our stars would have to be increased (or decreased in SC-stars) systematically by 0.3 to 
0.5 dex depending on the individual element. In that case SC- and C(N)-stars s-element 
overabundances can be estimated at similar levels within the error bars. This fact is a clear 
example of the large uncertainties involved in the chemical analysis of AGB stars. 

As to the class of MS/S stars, which are almost as close to solar metallicity as the 
present sample of C-stars, the average <[ls/Fe]> is +0.7, with a typical uncertainty of 
0.25 dex (Busso et al. 2001). Within the respective uncertainties they have therefore 
enhancements very similar to those of C(N) stars in the present paper. With respect to 
the warmer Ball giants in the disk, which show <[Fe/H]> = —0.3, the subsample with 
[Fe/H] > —0.3, amounting to 8 stars in the selection made by Busso et al. (2001), gives 
<[ls/Fe]> = +0.85. Actually, the Ba stars are considered to be extrinsic AGBs, i.e. binary 
systems in which the present giant was enriched in s-elements by mass transfer from the 
primary AGB star. This means that the original AGB envelope must have been even more 
enriched in low-mass s-elements, by a factor of at least 0.2 dex. The present sample of 
C-stars as well as that of similar s-enriched stars is not an ample one and the average values 
should be taken with some caution. Furthermore, both the MS/S and the C-stars show 
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more complex spectra crowded by molecular lines than the warmer Ba II giants, which are 
also less C-rich. Summing up, we shall have to wait for more ample statistics and refined 
model atmospheres of the C-rich stars to make further inferences in these comparisons. 

Three stars in the sample, V CrB, VX Gem and lY Hya, deserve special attention. 
The first two stars are Mira variables with periods of 357 and 379 days, respectively. For 
the Mira variable V CrB, Kipper (1998) derived [Fe/H] = -2.12 and [Is/Fe] ^ +1.5 and 
concluded that the abundance pattern in this star is similar to that of other late-type 
CH-stars. The values we derived in this star for these figures are very different (see Table 2). 
We do not have an easy explanation for this discrepancy, as we used the same model 
atmosphere (kindly provided by T. Kipper). Because of T^s variations in Mira variables, we 
have analyzed this star using other Teg values as proposed in the literature (see discussion in 
Kipper's paper). However, the best fit to the observed spectrum is obtained with the value 
adopted here, Tgs — 2250 K. Inspection of Figure 3 in Kipper (1998) makes clear that the 
heavy element abundances derived by spectral synthesis just fit the cores of the absorption 
lines without reproducing consistently the global behavior of the observed spectrum. We 
thus believe that the line list used in Kipper's analysis of V CrB does not contain important 
CN and C2 absorption, what unavoidably leads to the derivation of larger [X/Fe] ratios. 
This could be the reason for the important differences between the two analyses. 

Also for the Mira variable VX Gem the spectroscopic abundance analysis is rather 
complex. This star shows line-doubling in the stronger fines (Na I D doublet, 7698 A K I, 
6708 A Li I). The hydrogen lines (q;,/5,7) as well as the 5980 A Mg I line are in emission. 
The Ca II infrared triplet shows P-Cygni profiles. All these features indicate the existence 
of a circumstellar envelope. On the other hand, T^s for this star is rather uncertain. We 
have adopted 2500 K according to the analysis of Grocncwegen et al. (1989), but Ohnaka 
& Tsuji (1996) derived 3050 K from infrared photometry and classified this star as a 
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SC-star. Indeed, effective temperatures in Mira variables change during the pulsational 
phase. However, we have found a better agreement between the theoretical and the observed 
spectrum in all the spectral ranges analyzed when adopting T^s = 2500 K. The abundances 
derived in this star therefore have to be considered with caution. In fact, although the Rb I 
line is clearly visible in its spectrum we prefer not to estimate the [Rb/M] ratio because of 
the probable contamination with a circumstellar Rb line. 

The situation concerning lY Hya is even more confusing. This star is a very red object 
with a high mass-loss rate, 4 x 10^® Mq yr~^ (Jura & Klcinmann 1989). Furthermore, 
lY Hya is one of the few Galactic super Li-rich C(N)-stars, with log (Li/H) -|- 12 ~ 4.0 
(Abia, Pavlenko, & de Laverny 1999). The atmosphere parameters for this star are quite 
uncertain since no direct determination exists in the literature. Evidence of this is the fact 
that we do not find a unique CNO abundance choice (within error bars) that would give 
a good fit for all the spectral ranges analyzed here. Thus, because of the very uncertain 
abundance ratios derived in V CrB, VX Gem and lY Hya, we have plotted their abundance 
ratios with smaller symbols in the figures below, to indicate that they have a lower weight 
in our analysis. 

4.1. Technetium 

As an additional hint to the abundance pattern found in our stars we searched 
for the presence of technetium. The search was based mainly on the analysis of the 
intercombination Tc I line at 5924.47 A (see Paper 1 for a detailed explanation of the 
procedure followed in the analysis of this Tc blend) . Unfortunately, the stronger resonance 
Tc fines around 4260 A are inaccessible in most C-stars because of the strong fiux depression 
in these stars below ~ 4400 A. In the last column of Table 2 we have indicated whether or 
not Tc is identified for each star. A "yes" entry means that the best fit to the 5924 A blend 
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is obtained when a Tc abundance different from zero is used in the spectral synthesis^. 
With a "no" entry we mean that a synthetic spectrum with no-Tc does not significantly 
differ from another one computed with a very small Tc abundance. In some stars the best 
fit to the Tc blend is obtained using a small Tc abundance but not strictly zero. These stars 
have been quoted with a "doubtful" (dbfi) entry since in our judgment the evidence favours 
the absence of Tc. We have to note however, that a "no" in Table 2 does not necessarily 
mean that Tc is absent since, as the 5924 A Tc I line is much weaker and in a more crowded 
zone than the resonance lines near 4260 A, a detailed study of the latter might eventually 
reveal the presence of Tc. Nevertheless, in two stars for which it was possible to analyze 
both the 5924 A and 4260 A spectral ranges (VX Gem and W Ori), we reached the same 
conclusion for both hues (see Table 2). If we exclude the stars quoted as doubtful, 50% of 
the stars in Table 2 show Tc, the others are compatible with a non-detection. Interestingly 
enough, the same figure is found in MS/S stars with/without Tc (Smith & Lambert 1985, 
1986, 1990). If one applies the same rules currently accepted for the MS/S giants to our 
C-rich objects, then Tc-no C(N)-stars would be extrinsic AGBs (candidates to have a WD 
companion), with extra dilution of s-element abundances related to mass transfer from a 
companion. However, we do not find any correlation between the level of heavy-element 
enhancements and the presence of Tc. In fact, for the Tc-yes and Tc-no stars we found 
<[ls/Fe]> = +0.53 and <[ls/M]> = -^0.48, respectively. 

As far as we know, there is no evidence of binarity for any of the stars studied here; 
despite this, we can take a further step in order to identify the extrinsic (binary) or 
intrinsic nature of our C-stars by studying their IR colours in the same way as Jorissen 
et al. (1993) did for S-stars. Infrared colours provide a measure of the optical properties 

''Examples of the quality of the fits to the Rb and Tc blends can be appreciated in the 
corresponding figures shown in Papers I and II. 
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of circumstellar dust. Such dust is expected to be around intrinsic (non-binary) C-stars, 
and thus the presence or absence of dust emission might enable us to discriminate between 
intrinsic and extrinsic C-stars. Following the results of Jorissen et al. (1993) most Tc-rich 
S stars have IR excesses {R > 0.1), which would indicate that they owe their chemical 
peculiarities to intrinsic nucleosynthesis processes. On the contrary, Tc-deficient (binary) 
S stars usually show 7? < 0.1. To estimate this IR excess due to dust, the fiux ratios 
R — F(12/xm)/F(2.2//m) were computed. The 2.2 /xm infrared flux was derived from the K 
magnitude taken from the Two Micron Sky Survey (Neugebauer & Leighton 1969; Claussen 
et al. 1987) and the 12 /xm fluxes from the IRAS point source catalogue. The small sample 
of C-stars and the fact that our Tc analysis is based on a weak line in a crowded zone 
prevents us from reaching any flrm conclusion. However, we flnd that 72% of Tc-no stars in 
Table 2 have R < 0.1, while this flgure decreases to 45% in the Tc-yes stars. This may be 
an indication of some correlation in the same direction as found in S stars. A similar study 
in a larger sample of C-stars is clearly necessary. 

5. The stelleir and nucleosynthesis reference models 

In order to compare our observed abundances of C(N)-stars with predictions from 
current models of AGB evolution and nucleosynthesis, we use the s-process nucleosynthesis 
results from Gallino et al. (1998), Travaglio et al. (2001) and Busso et al. (2001), which 
are based on a large set of AGB stellar models, for stars of initial mass between 1.5 and 7 
Mq, made with the FRANEC evolutionary code (Straniero et al. 1995, 1997, 2001). The 
s-process nucleosynthesis calculations were performed with a post-process code and were 
extended over a wide range of metallicities. Specifically, we make use here of the models 
for 1.5 and 5 Mq stars, in order to show how different predictions about abundances and 
abundance ratios at the Sr-Y-Zr peak that can be obtained when the initial mass is assumed 
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to belong to the LMS or to the IMS range. Reimer's (1975) parameterization for mass loss 
was adopted (77 = 0.7 and 10, for the 1.5 Mq and the 5 Mq cases, respectively). When 
models are run at nearly solar metallicity, the one relevant for most of our sample stars, 
third dredge-up is found to start after a few thermal pulses. The parameter A (defined as 
AM(TDU)/AMh, where AM(TDU) is the mass dredged-up in one mixing episode and 
AMh is the extension in mass by which the H shell advances in an interpulse period) is on 
average 0.24 and 0.38 in the two cases. TDU is followed in detail pulse after pulse, mixing 
with the envelope material from the He intershell rich in ^^C and s-elements. In this way 
we can follow the stellar envelope composition evolving in time while the star ascends the 
TP-AGB phase. 

The activation of the ^^C neutron source requires some protons from the envelope to 
penetrate into the He intershell when TDU is established. Subsequently, proton captures on 
the abundant ^^C allow the formation of a small ^^C-pocket in the uppermost layers of the 
He intershell. All ^^C nuclei are then consumed by a-captures during the interpulse phase, 
while the intershell zone is in radiative conditions (Straniero et al. 1997). Successful models 
of the ^^C-pocket formation have been presented in the recent years, as a consequence of 
the application of diffusive or hydrodynamical simulations of the partial mixing occurring 
at the H-C interface (Iben k Renzini 1982; HoUowell k Iben 1988; Herwig et al. 1997; 
Cristallo et al. 2001). A general consensus is emerging that such partial mixing events do 
indeed occur: they might be stimulated by rotational shear (Langer et al. 1999). However, 
the above models are still very speculative and cannot be used to make general predictions 
suitable to perform a comparison with observations. Therefore, the amount of ^^C in the 
pocket and the distribution within it must still be considered as free parameters. Following 
the approach by Busso et al. (2001), we chose here to adopt a wide range of ^^C production 
in the pocket. We start from the reference choices made in Gallino et al. (1998) for LMS, 
and in Travaglio et al. (2001) for IMS: both choices were called standard (ST) in the 
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original papers. They correspond to a total ^^C mass of 4 x 10~^ Mq in LMS and of 4 x 
10"''' Mq in IMS. The decrease by an order of magnitude of the strength of the -"^^C neutron 
source in the pocket corresponds to the overall decrease by the same order of the mass 
of the He intershell, the TDU efficiency, and the interpulse period. Then we scale such 
abundances upward and downward, obtaining cases with efficiencies from a factor of two 
higher (STx2) to a factor of 12 lower (ST/12). 

The envelope composition is progressively enriched by TDU in C and s-elements, 
reaching the condition for a C(N) star (C/0 > 1) in the last TPs. For AGB stars of initial 
mass 1.5 Mq this occurs for any metallicity up to [Fe/H] — +0.2. Table 5 shows the 
predicted C/0 and ^^C/^^C ratios in the envelope after the last TDU episodes where the 
condition C/0 > 1 is met for AGB models of 1.5, 3 and 5 Mq and for the three metallicities 
[Fe/H] = +0.2, 0.0 and —0.3. For each mass, in the third column we also report the 
remaining envelope mass at the given TP. With the adopted mass loss prescriptions, at 
[Fe/H] = +0.2 we do not find C/0 > 1 in the envelope of AGB stars of 3 and 5 Mq, 
although we are not far from reaching this condition. Our sample stars span a range —0.3 
to +0.3 in [Fe/H], but the estimated uncertainty of the metallicity and, even more, the 
observational uncertainty in the mass loss rate does not allow us to make a more detailed 
analysis such as, for example, to exclude on that all stars with nominal [Fe/H] higher than 
solar are of initial mass below 3 Mq. 

6. Discussion 

6.1. The neutron density and the mass of C stars 

Figures 1 and 2 show the run of observed abundance ratios [X/M], for Rb, Sr, Y and 
Zr, compared with AGB model predictions for stars of initial mass 1.5 and 5 Mq and 
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different metallicities. The plotted curves pick up the envelope composition when the 
photospheric C/0 ratio reaches unity, because we must fulfill the general requirement 
that our stars are carbon rich with C/0 ratios in the range 1.01 to 1.20. The s-element 
abundance predictions for AGB stars of M = 3 Mq are almost identical to the M — 1.5 
Mq case. An exception is [Rb/Fe], which is 0.1 dex higher for all choices of the ^^C-pocket. 
Different mass loss prescriptions or details on TDU pulse after pulse would not affect much 
the results, provided C/0 = 1 is reached in the envelope. 

Table 6 gives the Rb, Sr, Y, and Zr expectations at C/0 =1 for the various AGB 
models and different ^^C-pocket prescriptions. All sample stars are reproduced by our 
theoretical models in the range ST/1.5 to ST/3. Figures 1 and 2 show only a weak tendency 
of the curves to be non-monotonic versus [M/H] in the range of the bona fide observed 
C(N) stars. For the 5 Mq model a decrease in s-element abundances with decreasing 
metallicity appears. This predicted trend of envelope abundances is quite different from 
what has recently been shown in Busso et al. (2001) for the same AGB stellar models, 
where in particular a steep increase of [Is/Fe] resulted with decreasing [Fe/H], up to a 
maximum at around +1.2, (which is reached at [Fe/H] ~ —0.5 for the ST case), followed 
by a decrease at lower metallicities. The discrepancy is only apparent and deals with the 
particular constraint to be meet here, i.e. with the mentioned fact that in our sample 
C(N) stars, C/0 is close to 1 (see Table 1). The trends shown in Busso et al. (2001) are 
instead representative of the last occurrence of TDU in the models. For a 1.5 Mq star, in 
our calculations this corresponds to the 17th TDU episode; for solar metallicity this implies 
C/0 = 1.44. 

In AGB stars of lower metallicity, given that the ^^C production in the top layer of 
the He intershell is always close to 0.2 in mass fraction, the carbon enrichment in the 
envelope is more easily established. Therefore the envelope C/0 ratio increases steadily 
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with decreasing metallicity, and is moderated only by the fact that at low [Fe/H] the stars 
condense from an 0-enhanced interstellar medium (Wheeler, Sneden & Truran 1989). At 
[Fe/H] = —1.3 (close to the estimated metallicity of V CrB) in a 1.5 Mq model after 17 
TDUs, C/0 = 11.4 and ^^C/^^C = 2221. For such a metallicity, a C/0 ratio close to unity 
is found at the very early occurrence of TDU (with a corresponding ^^C/^^C = 135), when 
the s-process abundances have been barely affected by neutron captures in the He intershell. 
We underline that this finding should be considered as a general rule for C(N) stars, as 
observations always reveal a moderate C enrichment. This may imply some selection effect: 
intrinsic AGB C-stars with higher C/0 ratios should exist, as high values of C/0 have been 
observed in both post- AGB stars reaching the pre-planetary nebula stage (Reddy, Bakker, 
& Hrivnak 1999) and in extrinsic C-rich objects originated by mass transfer phenomena 
(Vanture 1992c). However, an excess of carbon is immediately translated into a copious 
production of carbon-rich dust and into a high opacity of the circumstellar environment, to 
the point of hiding the photosphere at optical wavelengths. This is the case of some infrared 
carbon stars observed by IRAS (van Loon et al. 1998; Marengo et al. 1999). Therefore, 
unhke MS-S giants or the extrinsic Ba stars, which display different levels of C enrichment, 
the visible C(N) stars represent a sort of snapshot taken at a particular moment during the 
TP- AGB evolution, always representing the first appearance of a carbon-rich composition. 

Figures 1 and 2 show, however, that non-negligible differences in the s-process 
predictions at C/0 —1 from a LMS and from an IMS exist. The IMS models show a 
considerably narrower strip of predicted abundances when the amount of ^^C burnt in the 
pocket is varied as compared with the LMS case. This is true especially for [Rb/Fe]. These 
IMS predictions, as a whole, seem scarcely compatible with most sample stars. 

In Figures 3 and 4 we plot the elemental ratios [Rb/Sr], [Rb/Y] and [Rb/Zr] versus 
metallicity. Of the three panels in the figures, the one yielding the clearest distinction 
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is the second one, [Rb/Y]. According to the discussion of § 2, the fact that IMS models 
(where the ^^Ne source is active) compare less favourably with the data is an indication 
that most neutron fluxes at the origin of the observed s-element enhancements come from 
the ^^C source, and consequently that most C(N) stars are of low initial mass. 

Eventually, in Figure 5 we use [Is/M] as the average abundance of [Y/M] and [Zr/M], 
after excluding Sr, whose abundance is more uncertain. The entire sample of stars having 
a measurement for the three elements: Rb, Y and Zr, are (with the exception of lY Hya) 
compatible with a metallicity [M/H] = 0.0 ± 0.30 and actually fulfills this requirement 
within 1 (7. Figure 5 considers only model sequences above [Fe/H] = —0.3. Along the 
theoretical curves, the predicted [Rb/Y] increases with increasing [Fe/H]; exception is the 
ST case for M = 1.5 M©, for which the lower [Fe/H] corresponds to the lower [Is/Fe]. 
The plot immediately displays the characteristics discussed above: LMS models occupy 
the region of negative [Rb/Y], IMS models stay on the opposite side, and the partial 
overlapping of figures 1 and 2 is avoided. 



6.2. Carbon isotopes and related problems 

Figure 6 shows in graphic form the C/0 observed ratio versus ^'^C/^^C. Current 
stellar models yield a small reduction in C/0 (by 30 %) and a consistent reduction of the 
i2q^i3q ratio (by a factor 4) during the first ascent of the red giant branch, caused by 
the first dredge up mechanism (FDU). The reduction by a factor of 4 from the solar ratio 
(^^C/^'^C = 89) is obtained after decreasing by one third of the initial ^^C abundance and 
through an increase by a factor ~2.5 of ^^C, both changes being due to proton captures 
on ^^C during the main sequence in the deep regions of the radiative H-rich envelope (e.g.. 
Dearborn 1992). Note that the resulting ratio after the FDU is independent of the initial 
i2q^i3q value. The same FDU results are obtained for different initial metallicities, even 
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when one accounts for the fact that, in Galactic evolution, A'"(^^C) scales as N{Fe), while 
N(^^C) scales as N{Fey (Iben 1977; Boothroyd & Sackmann 1988). 

However, at least for masses smaller than about 2 Mq the canonical reduction of the 
i2q^i3q ratio at FDU down to 20 — 25 is insufficient to explain the observations in red 
giants (Gilroy & Brown 1991; Charbonnel, Brown, & Wallerstein 1998; Gratton et al. 2000). 
A similar discrepancy with model predictions is shown by oxygen and carbon isotopes 
measured in dust grains of clear AGB circumstellar origin, recovered from meteorites. 
This is true both for C-rich (Anders & Zinner 1993) and for 0-rich grains (Nittler et 
al. 1994; Hutcheon et al. 1994). Explanation of this evidence involves a relatively slow 
process of mixing occurring during the RGB in the radiative region located between the 
bottom of the convective envelope and the H- burning shell, perhaps induced by rotation 
(Sweigart & Mengel 1979). This process is sometimes called Cool Bottom Process, or 
GBP (Wasserburg, Boothroyd, & Sackmann 1995). GBP has been shown to yield quite 
naturally enhancements of ^^G and of ^^O and reductions in ^^O. Several phenomenological 
models have been presented for such phenomena (Boothroyd, Sackmann, & Wasserburg 
1994; Denissenkov & Weiss 1996; Weiss, Denissenkov, & Gharbonnel 2000). The general 
agreement with observations has led to acceptance that such circulations should indeed 
be active, thus further reducing in LMS red giants the ^'^C/^^C ratio established after the 
FDU. The exact extent of the depletion as a function of mass for each metallicity is, to date, 
rather uncertain, as the requirement of sampling low mass objects has often led to obtain 
the measurements in metal poor clusters, so that the two parameters (mass and metal 
content) are coupled and difficult to separate. Recently, Charbonnel, Brown, & Wallerstein 
(1998) have started a project on field stars whose parallax was accurately measured by 
HIPPARGOS, and a clearer view is now emerging, in which giant stars of the old Galactic 
disc population (at metallicities [Fe/H] ~ —0.6) can reach -'^^G/-'^^G ratios close to 7. For 
more metal-rich stars a reasonable estimate might be 10 — 15 (Gilroy & Brown 1991). Not 
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having a suitable model for CBP, we chose this estimate as the ^^C/^^C ratio present in 
early AGB phases of LMS before the onset of thermal pulses: the adopted tentative value 
is 12. For stars more massive that ~ 2.5 Mq, we assumed ^^C/^^C = 24. 

When we then add ^^C to the envelope by TDUs along the AGB, the isotopic ratio 
raises. At C/0~ 1 a pre- AGB value of around 25 would translate into a final value a 
bit higher than the solar isotopic ratio (^^C/^^C=89). Indeed, after a few TPs, when 
TDU is operating, a nearly constant mass fraction of ^^C , X(^^C) ~ 0.23, is built in the 
He intershell by partial He burning in the TP. This is independent of the core mass, i.e. 
on the AGB initial mass. About identical results are obtained by different authors (see 
e.g., Iben 1976, 1982; Boothroyd & Sackmann 1988; Straniero et al. 1995, 1997; Mowlavi 
1999). Only a marginal production of ^^O occurs in the TP, so that the oxygen abundance 
in the envelope is not affected by TDU. Also the ^^C abundance in the envelope is not 
affected by TDU, since no ^^C nuclei are left in the He intershell, while the mass of the 
He intershell cumulatively mixed with the surface is a small fraction of the envelope. 
Under those conditions, it is easy to evaluate the ^^C/^^C in the envelope at C/0 = 1. 
It results: ^^C/^^C ~ 3.6 x (^^C/^^C )fdu- This relation is independent of metallicity if 
both ^^C and ^^O are assumed to scale as Fe. Actually, for [Fe/H] = —0.3 we introduced 
a modest a-enhancement in the initial ^^O abundance with respect to a solar-scaled 
metallicity composition, according to an average observational trend [0/Fe] = —0.5 [Fe/H] 
(cf. Wheeler, Sneden, & Truran 1989). This implies an increase by about 30% of the 
predicted ^^C/^^C ratio at C/0 = 1. 

As shown in figure 6, only two stars, V Aql and W Ori, show the expected ratio at C/0 
~ 1, without requiring any CBP to operate. Those two stars should have a mass higher 
than 2.5 Mq, perhaps around 3 Mq. When we consider the majority of our sample stars, 
which have previously been shown to be of low mass, we immediately see that CBP must 
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have reduced the ^^C/^^C ratio. Actually, even a CBP operating only on the RGB would 
be insufficient, as according to the previously given rule, starting from (^^C/^^C)fdu ~ 12 
we would have final values of 43 (see Figure 6). Hence CBP should be active also later, 
in the early AGB phase, to operate changes larger than we assumed and thus explain the 
later part of the ^^C/^^C distribution. The possibility exists that CBP operates also in the 
TP- AGB phase, though calculations have not been performed yet. A different efficiency of 
rotationally induced meridional mixing (Sweigart & Mengel 1979) in different stars might 
also accomplish this job. 

Furthermore, in the above discussion we assumed that all sample stars are intrinsic 
AGBs. This actually may not be true. A considerable fraction of C(N) stars might be 
extrinsic AGBs belonging to binary systems (see § 4.1). In this case, the observed star could 
be ascending the early AGB phase, having being enriched in carbon and s-process elements 
by mass transfer from the more massive companion while this last was on the AGB. In this 
case FDU on the companion star can further reduce the ^^C/^^C ratio ^. 

7. Conclusions 

In this paper we have presented new high-resolution spectroscopic measurements 
for carbon- rich AGB stars belonging to the C(N) type. We have focused our attention 
on s-elements around the ^^Kr reaction branching of the s-process path, showing that 
the comparative analysis of Rb, Sr, Y, and Zr yields information on the neutron density 
prevailing during the neutron capture nucleosynthesis processes occurring in the He 

^Note, that actually there is a controversy on the typical ^^C/^'^C ratio in C(N)-stars (see 
Lambert et al. 1986; Ohnaka & Tsuji 1996-98; de Laverny & Gustafsson 1998-99; Schoier & 
Olofsson 2000). 
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intershell. Despite the large uncertainty existing in the derivation of element abundances 
in these stars, comparison with AGB nucleosynthesis models, which were already capable 
of explaining the abundances in other intrinsic and extrinsic classes of AGB stars, confirms 
a general agreement between theoretically predicted photospheric compositions and 
observations. The main implication of such an analysis is that most C(N) stars are of low 
mass, experiencing s-process nucleosynthesis phenomena dominated by the neutron source 
provided by a captures on -"^^C in radiative intershell layers. When the ^^C/^^C ratios are 
considered, their generally low values require the activation of slow circulation mechanisms 
below the formal convective border of the envelope, bringing material close to the H-burning 
shell and therefore increasing the abundance of ^^C. These mechanisms are indeed known to 
occur in stars below about 2 Mq. this implies, therefore, the conclusion that the majority 
of C(N) stars are of small mass. Exceptions might exist in a very few cases, though this 
speculation is highly uncertain and should be verified by more extensive measurements of 
s-process elements, especially at the Ba-peak. 
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Tabic 1. Data for Program Stars 



Star Spec. lype Teff Ref.^ C/0 ^'^G/^^C Ref.^ 



AQ And 


C5,4 


2970* 


2 


1.02 


30 


1 


AWCyg 


C3,6 


2760 


2 


1.03 


21 


2 


EL Aur 


C5,4 


3000 


6 


1.07 


50 


2 


HK Lyr 


C6,4 


2866 


2 


1.02 


10 


2 


IRC -10397 


N 


2600 


6 


1.01 


20 


2 


lY Hya 


N 


2500 


6 


1.02 


15 


2 


LQ Cyg 


C4,5 


2620 


2 


1.10 


40 


2 


RX Set 


C5,2 


3250 


2 


1.04 


60 


2 


S Set 


C6,4 


2895 


1 


1.07 


45 


3 


SY Per 


C6,4 


3070 


3 


1.02 


43 


2 


SZ Sgr 


C7,3 


2480 


2 


1.03 


8 


2 


TY Oph 


C5,5 


2780* 


2 


1.05 


45 


2 


U Cam 


C3,9 


2670 


2 


1.02 


40 


2 


UU Aur 


C5,3 


2825 


1 


1.06 


50 


3 


UV Aql 


C6,2 


2750 


2 


1.005 


19 


2 


UX Dra 


C7,3 


2900 


1 


1.05 


26 


2 


V Aql 


C6,4 


2610 


1 


1.16 


90 


2 


V CrB 


C6,2 


2250 


4 


1.10 


10 


4 


V Oph 


C5,2 


2880* 


2 


1.05 


11 


1 


VX Gem 


C7,2 


2500 


5 


1.01 


60 


2 
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Table 1 — Continued 



Star Spec. lype T^s Ref.^ C/0 ^^C/^^C Ref.^ 



V460 Cyg 


C6,3 


2845 


1 


1.06 


61 


4 


V781 Sgr 


N 


3160 


2 


1.02 


35 


3 


W Ori 


C5,4 


2680 


1 


1.20 


79 


4 


WZ Cas 


C9,2J 


3140 


3 


1.005 


4 


3 


Z Psc 


C7,2 


2870 


1 


1.01 


55 


4 



^References for T^s- (1) Lambert et al. (1986); (2) Ohnaka & Tsuji (1996); 
(3) Dyck, van Belle & Benson (1996); (4) Kipper (1998); (5) Groenewegen et 
al. (1998); (6) For these stars we have estimated the effective temperature by 
comparing their spectra with the temperature sequence spectral atlas for C-stars 
created by Barnbaum et al. (1996). 

''Source for the ^^C/^^C ratio: (1) Ohnaka & Tsuji (1996); (2) Derived in this 
work; (3) Lambert et al. (1986); (4) Kipper (1998) 



-44- 



Table 2. Abundances Derived in Program Stars 



Star 


[M/H] 


[Rb/M] 


[Sr/M] 


[Y/M] 


[Zr/M] 


Tc 


AO And 




u.o 


n 6 


n 7 


u.o 


llU 


AW uyg 


n n 


u.z 


n A 


u.o 


u.o 


yes 


i-jij -TiUI 


— u.uo 


n /I 


u.o 


n 7 

U. ( 


n 8 
u.o 


no 


WT^ T -iry 

riiv i-iyr 


n 1 n 


n 
u.z 


u.o 


n 7 

U. / 


n 7 
u. / 


yes 


lil*^ — iUoy/ 


n 1 n 

U.iU 


n 1 

U.i 


u.o 








TV Ha^ci 


— U.oU 


u.o 


u.o 


u.o 


1 n 

i.U 






u.zo 


u.z 


u.o 


u.o 


n 7 

U. / 


no 


rlA. OCX 


— u.uo 


n /I 

U.4 


n /I 
u.4 


u. / 


n 8 
u.o 


no 


ij OCX 


U.Ul 


u.o 


u.o 


n 7 
u. / 


u.o 


yes 


Pot 

o 1 r er 


— U.oU 


u.i 


u.o 






QOUDI 


oZj Ogl 


— U.U^ 


n 1 

u.i 


n A 

u.4 


n 8 
u.o 


n 8 


QOUUi 


TV r»T-vVi 
i I wpn 


U.iU 


u.z 


U.D 


n 7 

U. / 


U.D 


doubf 


TT PciTTl 


n riQ 
— u.uy 






u.o 


n A 
u.4 


no 


TTTT Anr 


u.uu 


n A 
u.4 


u.4 


u.o 


u.o 


no 


UV Aql 


-0.06 


0.6 


0.5 


0.7 


0.9 


yes 


UX Dra 


-0.20 


0.0 


0.8 


0.8 


0.8 


yes^ 


V Aql 


-0.05 


0.3 


0.5 






no 


V CrB 


-1.35 


-0.2 


0.3 


0.4 




no^ 


V Oph 


0.0 


-0.2 


0.2 


0.3 


0.4 




VX Gem 


-0.15 




0.6 


0.5 


0.6 


yes^ 


V460 Cyg 


-0.04 


0.4 


0.2 


0.4 


0.4 


yes 
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Table 2 — Continued 



Star 


[M/H] 


[Rb/M] 


[Sr/M] 


[Y/M] 


[Zr/M] 


Tc 


V781 Sgr 


0.10 


0.1 


0.5 


0.3 


0.2 


yes 


WOri 


0.05 




0.0 


0.3 


0.1 


no^ 


Z Psc 


-0.01 


0.9 


0.9 


1.0 


0.8 


yes 



'^Little et al. (1987) classified SZ Sgr as possible star showing 
Tc and reported Tc detection in UX Dra. Kipper (1998) sets 
an upper limit to the Tc abundance in V CrB (log (Tc/H) +12 
< —0.5) from the analysis of the 5924 A Tc 1 line. However, 
our analysis is compatible with no Tc. For VX Gem and W Ori 
we were able to get spectra in the 4260 A region whose analysis 
confirmed the figure found from the 5924 A Tc I analysis. 



Table 5. C/0 and ^^C/^^C Predictions in the Envelope 







i.O Mq 






Q T\/r 

6 Mq 






c; T\iT 

Mq 




Metallicity 






/i/f 






/i/f 

Menv 






^^env 


[Fe/Hl= +0.2 


1.00 


43 


0.28 
















1.05 


45 


0.24 














[Fe/H]= +0.0 


1.06 


46 


0.38 


1.00 


91 


1.05 


1.00 


92 


1.23 




1.17 


51 


0.34 


1.04 


94 


0.98 


1.04 


95 


1.16 




1.27 


55 


0.31 


1.08 


98 


0.91 


1.08 


99 


1.09 




1.37 


60 


0.28 


1.11 


101 


0.83 


1.12 


103 


1.02 




1.45 


63 


0.24 


1.14 


104 


0.76 


1.16 


107 


0.95 
















1.21 


111 


0.88 
















1.26 


116 


0.80 
















1.32 


122 


0.73 


[Fe/H]= -0.3 


1.01 


62 


0.44 


1.01 


131 


1.31 


1.01 


132 


1.58 




1.16 


72 


0.41 


1.08 


140 


1.24 


1.05 


137 


1.51 




1.31 


81 


0.38 


1.13 


146 


1.18 


1.09 


143 


1.44 




1.47 


92 


0.34 


1.18 


154 


1.11 


1.13 


149 


1.37 




1.61 


101 


0.31 


1.23 


160 


1.05 


1.18 


155 


1.30 




1.76 


110 


0.28 


1.28 


167 


0.98 


1.23 


162 


1.23 




1.86 


117 


0.24 


1.34 


174 


0.91 


1.28 


169 


1.16 










1.38 


180 


0.83 


1.33 


176 


1.09 










1.44 


187 


0.76 


1.39 


184 


1.02 



-47- 



Table 6. Abundance Ratios at C/0= 1 for different AGB models and ^^C-pocket 

prescriptions 



Model i3Q_pQ(,]^g^ c/0 [Rb/Fe] [Sr/Fe] [Y/Fe] [Zr/Fe] 



1.5 [Fe/H]= +0.2 ST 1.00 43 0.51 0.91 0.82 0.70 

ST/1.5 1.00 43 0.31 0.61 0.52 0.39 

ST/3 1.00 43 0.09 0.16 0.10 0.05 

1.5 Mq [Fe/H]= +0.0 ST 1.06 46 0.56 1.06 0.96 0.82 

ST/1.5 1.06 46 0.37 0.75 0.65 0.54 

ST/3 1.06 46 0.12 0.28 0.20 0.12 

1.5 Mq [Fe/H]= -0.3 ST 1.07 46 0.39 1.12 1.12 1.08 

ST/1.5 1.07 46 0.37 0.95 0.86 0.73 

ST/3 1.07 46 0.16 0.46 0.37 0.29 

3 Mq [Fe/H]= +0.0 ST 1.00 91 0.71 1.05 0.98 0.85 

ST/1.5 1.00 91 0.51 0.75 0.68 0.57 

ST/3 1.00 91 0.20 0.29 0.22 0.14 

3 Mq [Fe/H]= -0.3 ST 1.07 97 0.47 1.14 1.15 1.10 

ST/1.5 1.07 97 0.46 0.96 0.89 0.76 

ST/3 1.07 97 0.22 0.48 0.40 0.32 

5 Mq [Fe/H]= +0.0 ST 1.00 92 1.33 1.11 1.08 1.29 

ST/1.5 1.00 92 1.24 0.93 0.87 1.02 

ST/3 1.00 92 1.04 0.71 0.63 0.73 

5 Mq [Fe/H]= -0.3 ST 1.01 93 0.92 0.91 0.96 1.29 

ST/1.5 1.01 93 1.02 0.89 0.89 1.12 

ST/3 1.01 93 0.96 0.67 0.61 0.74 
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Fig. 1. — Derived abundance ratios [X/M] vs. [M/H] in our sample of C-stars (dots) com- 
pared with the corresponding theoretical predictions (hues) for a 1.5 Mq TP-AGB star when 
the envelope reaches C/0 > 1 for different parameterizations of the ^^C pocket (see text). A 
typical error bar is shown. In the stars where a very few number of lines are used the error 
bar is certainly larger. Note, furthermore, that the data points for lY Hya, VX Gem and V 
CrB are plotted with smaller symbols to indicate that the abundance ratios derived in them 
are very uncertain. 

Fig. 2. — As Figure 1 comparing with theoretical predictions for a 5 Mq TP-AGB star. 

Fig. 3. — Comparison of the observed [Rb/Sr,Y,Zr] ratios vs. [M/H] with different theoretical 
predictions for a 1.5 Mq TP-AGB star when the envelope reaches C/0 > 1. Note that again 
the data points for lY Hya, VX Gem and V CrB are plotted with smaller symbols. 

Fig. 4. — As Figure 3 comparing with theoretical predictions for a 5 Mq TP-AGB star. 

Fig. 5. — Comparison of the observed mean low-mass s-element (Y,Zr) enhancement (sig- 
nature of the neutron exposure) against [Rb/Y] (signature of the neutron density) with 
theoretical predictions for a 1.5 Mq (upper panel) and 5 Mq (lower panel) TP-AGB star for 
different parameterizations of the ^"^C pocket (lines). The theoretical predictions shown are 
for stellar models with C/0 ~ 1 and metallicity [Fe/H] > —0.3, which is the metallicity of 
most of the stars studied here. Only stars in the sample with metallicity above this value 
are plotted (see text). Note that several stars coincide in the same data point. 

Fig. 6. — C/0 vs. ^^C/^^C ratios derived in our stars. Note that most of the stars have a 
C/0 ratio only slightly higher than 1 and that a significant number show ^^C/^^C < 40. An 
average error bar is shown. 
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Table 3. Spectroscopic Data and Equivalent Withds 









AQ 


AW 


EL 


HK 


IRC lY LQ 


RX 


s 


SY 


sz 


TY 


X['^^ ) 




And 




Aur 


Lyr 


— 10397 Hva Cvp- 






Per 


Sgr 


Onh 

W JJii 


Cal 6798.48 


2.70 


-2.320 


200 




<260 


200 


150 


182 


126 


220 


240 




Til 4503.75 


2.13 


-0.840 






















Til 4759.63 


2.17 


-1.514 


















112 




Til 4801.94 


0.83 


-3.111 






208 


160 


<190 


160 


170 






160 


Til 4805.41 


2.34 


0.150 








260 


115 


162 


153 




145 


135 


Til 4812.89 


0.84 


-3.450 






148 




128 


120 










Til 4821.27 


2.16 


-0.940 






















Til 4870.12 


2.24 


0.518 


















186 




111 4yZD.14 


U.oi 


1 7n 




1 07 


070 






<zuu 










Til 4937.72 


0.81 


-2.254 






















Til 5238.53'" 


2.09 


-1.102 


221 




237 






187 








180 


Til 5238.58'" 


0.85 


-1.850 






















Til 5716.48 


2.29 


-0.700 


250 


200 








250 


227 


<300 


<230 


<240 


Til 7065.05 


1.46 


-2.609 






















V I 4748.50 


2.04 


-0.620 






















V I 5668.36 


1.08 


-1.030 




<250 


<300 


234 


185 ... 260 


<300 








<260 


V I 5670.85 


1.08 


-0.420 


277 






225 
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Table 3 — Continued 



AQ AW EL HK IRC lY LQ RX S SY SZ TY U 

A(A) x(eV) log 5/ And Cyg Aur Lyr -10397 Hya Cyg Set Set Per Sgr Oph Cam 



Fel 4445.48 


0.09 


-5.441 
















278 






Fel 4794.36 


2.42 


-4.050 




118 








92 










Fel 4809.93 


3.57 


-2.720 












100 


90 


40 






Fel 4813.11 


3.27 


-2.890 








77 


80 






88 




97 


Fel 4817.77 


2.22 


-3.530 


80 




140 


125 


110 


110 


119 


113 


100 


158 


Pel 4911.77 


3.92 


-1.687 
















51 






Fel 5848.12 


4.60 


-0.903 






















Nil 4756.51 


3.48 


-0.340 








120 








100 






Col 4813.50 


3.21 


0.050 




72 


89 


70 


70 


42 


51 


29 


47 




Rbl 7800.23 


0.00 


0.130 






















Sri 4811.88 


1.84 


0.070 








122 








150 ... 






Sri 4872.45 


1.79 


-0.200 






















Sri 7070.10 


1.84 


-0.150 






















Y I 4780.18 


1.88 


-0.910 




















<260 


Y I 4781.02 


1.39 


-0.460 


120 






150 








... 126 






Y I 4819.64 


1.36 


-0.480 


140 


100 


120 


140 


150 


145 


120 


... 155 


134 


140 


Y I 4879.65 


1.90 


-0.400 


















100 
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Table 3 — Continued 



A(A) 


AQ 


AW 


EL HK 


IRC lY 


LQ 


RX S SY SZ TY U 


x(eV) log gf And 


Cyg 


Aur Lyr 


-10397 Hya 


Cyg 


Set Set Per Sgr Oph Car 



Y I 8800.50 


0.00 


-2.240 




362 








Zrl 4772.32 


0.62 


0.040 








280 <270 




Zrl 4805.88 


0.63 


-0.420 


190 


220 210 


200 


200 185 


198 


Zrl 4815.03 


0.65 


-0.530 












Zrl 4815.62 


0.60 


-0.630 












Zrl 4828.05 


0.62 


-0.640 


230 


195 270 


210 ... 


180 200 ... 250 


202 


Zrl 7849.39 


0.68 


-1.300 




350 








Zrl 8389.50 


0.60 


-1.760 




360 








Tel 5924.72 


0.00 


-2.440 













"^hfs: hyperfine strueture included; ss: speetrum synthesis used; s: solar gf derived; BG: Biemont et al. (1978); CB: Corliss 
(1962); FMW: Fuhr et al. (1988); HL: Hannaford et al. (1982); MFW: Martin, Puhr, & Wiese (1988); NBS: Wiese & Fuhr 

Komarovskij (1991); No entry: Kurucz CD ROM No. 23 or VALD data bases. 

''These two Ti I lines eannot be resolved. The equivalent width shown for the 5238.53 A line corresponds to the wl 
Obviously, we consider the contribution of the two lines in the abundance analysis. 
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Table 3. Continued 



A(A) 


uu 


uv 


ux 


V 


V 


V 


vx 


V46() 


V781 


w 


z 


Aur 


Aql 


Dra 


Aql 


CrB 


Oph 


Gem 


Cyg 


Sgr 


Ori 


Psc 


Cal 6798.48 








218 




167 


125 


153 


<250 






Til 4503.75 






















140 


Til 4759.63 






96 


















Til 4801.94 


200 














<190 


151 


220 




Til 4805.41 


156 


116 












150 




150 


138 


Til 4812.89 


146 


115 


108 






115 








110 




Til 4821.21 






















125 


Til 4870.12 






180 


















Til 4926.14 
















<230 








Til 4937.72 




















238 




Til 5238.53 


<270 
















180 


186 


<265 


Til 5238.58 
























Til 5716.48 




264 


157 




170 


<200 


123 


250 


210 


<270 




Til 7065.05 
























V I 4748.50 




















109 




V I 5668.36 


250 


280 


286 










<270 


<260 


<285 


<280 


V I 5670.85 












<240 










208 
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Table 3 — Continued 



UU UV UX V V V VX V460 V781 W Z 

Aur Aql Dra Aql CrB Oph Gem Cyg Sgr Ori Psc 



A(A) 



Fel 4445.48 
Fel 4794.36 
Fel 4809.93 
Fel 4813.11 
Fel 4817.77 
Fcl 4911.77 
Fel 5848.12 
Nil 4756.51 
Rbl 7800.23 
Col 4813.50 
Sri 4811.88 
Sri 4872.45 
Sri 7070.10 

Y I 4780.18 

Y I 4781.02 



300 ... 

93 

67 

117 124 118 



54 
130 



86 



52 36 



110 135 



90 140 



Y 1 4819.64 145 143 155 
Y I 4879.65 



44 



20 



65 
63 



35 
90 



50 



105 
145 



65 



49 
61 



108 

80 
121 
120 



95 



50 



91 



136 

70 
90 
125 



43 



107 126 120 91 



75 



20 



114 130 



120 138 120 155 140 
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Table 3 — Continued 



UU UV UX V V V VX V460 V781 W Z 

A(A) Aur Aql Dra Aql CrB Oph Gem Cyg Sgr Ori Psc 



Y I 8800.50 

Zrl 4772.32 <350 <270 ... <350 

Zrl 4805.88 200 ... 135 194 210 185 170 182 190 

Zrl 4815.03 

Zrl 4815.62 

Zrl 4828.05 235 203 200 198 ... 200 185 196 195 

Zrl 7849.49 

Zrl 8389.50 

Tcl 5924.72 
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Table 4. Dependence of Derived Abundances on Model Atmosphere Paramel 



log Wa/A ATeff = ±200 K = ±0.5 kms-^ AC/0 = ±0.05 A^^C/^''^C= ±10 ACNO/ 



[Rb/H] 


ss^ 


±0.28 


^0.01 


^0.05 


^0.05 


[Sr/H] 


ss'^ 


^0.13 


^0.10 


TO.02 


^0.05 


[Y/H] 


-4.52 


±0.25 


^0.20 


TO.02 


±0.01 


[Zr/H] 


-4.39 


±0.15 


^0.30 


^0.05 


±0.01 


[M/H] 


-4.77 


±0.20 


^0.10 


^0.02 


±0.01 



'^Uncertainty due to a typical ±20 niA error in the equivalent widths. When spectral synthesis is used, t 
error in the placement of the spectral continuum. 

^The abundances of Rb and Sr (in most of the stars) were derived using spectral synthesis (ss), see tex 



